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Abstract. - 

A rigorous and simple perturbative proof of Luttinger's theorem is sketched for Fermi liquids 
in two and three dimensions. It is proved that in the finite volume, the quasi-particle density 
is independent of the interaction strength. The thermodynamic limit is then controlled to all 
orders in perturbation theory. 



For more than the past 40 years Fermi liquid theory has played a major role in the descrip- 
tion of condensed matter systems at low temperature. One of the most striking characteristics 
of a Fermi liquid is probably the presence of a Fermi surface, revealed by a discontinuity in the 
occupation number as function of the quasi-momentum at zero temperature. Luttinger found 
in his seminal paper [1] that the volume enclosed in the Fermi surface is independent of the 
interaction strength, if the density is held fixed. Assuming the existence of a formal expansion 
in powers of the interaction strength, Luttinger proved that the higher order contributions to 
the volume vanish order by order in the perturbation expansion. 

From a mathematical point of view, the proof given by Luttinger is unsatisfactory. The 
manipulations of conditionally convergent integrals that lead to Luttinger's result have to 
be controlled rigorously. One could introduce an infrared cut-off, in which case Luttinger's 
argument is mathematically valid, and carefully remove the cut-off in a second step, using 
renormalization group ideas. 

Recently, an alternative topological proof of Luttinger's theorem was given by Oshikawa [2]. 
On a finite torus, he considered the change in the total momentum in a particular direction 
under the adiabatic insertion of a unit quantum flux $o = hc/e through the hole of the torus. 
Expressing the change in the total momentum in terms of the volume of the Fermi sea, he 
showed that the volume of the Fermi sea is given (up to multiplication by a constant) by the 
particle density plus an integer, corresponding to the number of completely filled bands. As 
Luttinger did, Oshikawa assumed the existence of the thermodynamic limit, and completed 
the proof of Luttinger's theorem. 

Instead of following the argument of Luttinger, we propose here a simpler proof of Lut- 
tinger's theorem. We first observe that the number operator commutes with the Hamiltonian. 
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Thus, the ground state of the system in a finite spatial volume is an eigenstate of the num- 
ber operator, whose eigenvalues are natural numbers. It follows that at zero temperature, 
if the ground state is non-degenerate, the expectation value of the number operator is an 
integer. On the other hand, we show that for generic dispersion relations (we consider a single 
band system, as for example a Hubbard system far from half- filling), the expectation value at 
zero temperature of physical observables, like for instance the number operator, are analytic 
functions of the interaction strength in the finite volume. Obviously, continuous variation of 
the interaction strength cannot change the zero temperature expectation value of the number 
operator. We deduce that, in the finite volume, the particle density is given by a power series 
in the coupling constant, in which only the zero order Taylor coefficient is non-vanishing. 
If the thermodynamic limit of the expectation value of any physical observable is carefully 
controlled, the above statement remains true in the infinite volume. 

In order to avoid unphysical infrared divergences in the thermodynamic limit, we need to 
renormalize the single particle dispersion relation. In our renormalization scheme, the "in- 
teracting" dispersion relation, defining the physical Fermi surface, appears in the propagator. 
We introduce a counter-term that holds the Fermi surface fixed, and obtain a renormalization 
map between the interacting dispersion relation, and the bare dispersion relation, determined 
by the single-particle Schrodingcr equation. We now have all the ingredients needed to prove 
Luttinger's theorem, which is obtained by fixing the chemical potential to its physical value, 
both in the free fermion approximation, and in the interacting system. Comparison of the 
two results for the volume of the Fermi sea yields Luttinger's theorem. 

In the present work, the thermodynamic limit is controlled order by order in powers of 
the interaction strength. In that sense, our result, as well as the proof of Luttinger [1] or 
Oshikawa [2] are perturbative, in contrast to the recent proof [3] of the existence of a 2D 
Fermi liquid at zero temperature, where the convergence of the power series in the interaction 
strength is established. Details of the proof of the existence of the thermodynamic limit in 
the sense of formal power series will be presented in a forthcoming paper by Praz [4]. 

A necessary assumption for the validity of Luttinger's theorem seems to be gauge invari- 
ance, or equivalently charge conservation. The argument of Luttinger can be related to Ward 
identities. Oshikawa explicitly makes use of gauge invariance and we require particle number 
conservation. 

We consider many fermions on the finite lattice of size L d with periodic boundary condi- 
tions, where the size of the unit cell is one. The system will be described by the many body 
Hamiltonian H = Ho + Ay, where the single-particle Hamiltonian is diagonal in momentum 
space, 

H = L- d J2^Ha^- (1) 

k,<T 

Here k is in the dual lattice, and a G {f, 1} is the spin. e(k) is the one-particle, single-band 
dispersion relation. Strictly speaking, the dispersion relation should depend on the size of 
the system. In order to lighten the notations, we work directly with the dispersion relation 
e(k), defined on the first Brillouin zone, neglecting unimportant terms of order l/L. Precise 
assumptions on the properties of e(k) are given below. The creation and annihilation operators 
c k<r ano - CkCT satisfy the fermionic anticommutation relations {c^ a ,c pT } = (2nL) d 5kp5 r7T . The 
pair interaction between fermions is given by 

v = L -3d ^ 6 kl +k 3 ,k 2 +k 4 ?VT(ki - k 3 )c+ CT c+ T c k3(T c k4T (2) 

ki,...,k 4 
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where we shall assume that the potential v(k) is local, i.e., is a smooth function of k. Our 
proof can trivially be extended to potentials depending smoothly on the frequency. 

The thermal expectation value at zero temperature of any polynomial in creation and 
annihilation operators A is 

tr p[ L }A 

= — (3) 

where P^ is the projection on the ground state of Hq — fiN + XV where the number operator 
N = L~ d c^Cka, and ^ is the chemical potential. In the thermodynamic limit L — > oo, 
the naive power expansion of (A)l,x in A does not exist. Precisely, the coefficients of X n in 
the power expansion diverge as L — > oo for n > 3. The renormalization procedure cures these 
(unphysical) infrared divergences order by order in A. 

In the present work, we are concerned with the perturbative theory only. We will not 
be interested in the convergence of the whole perturbative series, which is related to (the 
absence of) instability in the ground state. We thus consider only formal power series. /(A) = 
Sn>o fn^ n is a well-defined formal power series in A if for all n, |/„ < oo. 

Infrared divergences appear because the bare propagator used in the expansion has a 
pole on the bare Fermi surface, while the physical pole lies on the interacting Fermi surface. 
We propose a renormalization scheme which allows a diagrammatic expansion around the 
interacting Fermi surface. Formally, we split the bare dispersion relation in a renormalized 
dispersion relation E(k; fi, A), and a counter-term -ft"(k; A) such that e(k) — fj, = E(h, /i, A) + 
K(k; A). The counter-term is fixed by the renormalization condition, which guarantees that 
the interacting Fermi surface is determined by the zero set of E(h, fi, A). Unfortunately, the 
renormalized dispersion relation cannot be determined directly. We thus begin with a generic 
dispersion relation e(k, jS), and consider the model defined by the Hamiltonian 



H ( L ) =L -dJ2e^ (k)c+ ff c kCT + XV + 3C< L > , (4) 

where 



k,cr 



[ e(k), if |e(k)| > 4tt/L, 

e^(k) = \ sgn (e(k))f , if |e(k)| < 4n/L, (5) 
[ 47r/L, if e(k) = 

is a regularized dispersion relation (the /x dependency has been dropped for simplicity), and 
%( L ) = L~ d a K^ L \k; e( • , fi), X)c^ a Ck a is the counter-term, to be determined by the (dis- 
cretized) renormalization condition. We then control the thermodynamic limit of the expec- 
tation value of an operator A, defined by 



tr yin 



(L) 
A, i 



E^(A;e,X) = —-^ (6) 

where 11^ is the projection onto the state of of minimum energy of the Hamiltonian Q. In 
the infinite volume, the occupation number of this model has a jump on the surface defined 

by 

£ e , M = {k: e(k,/i)=0}, (7) 
which is identified as the interacting Fermi surface. If 



e(k) — M = e(k, ^t) + K (k; A, e( • , fi)), 



(8) 
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we can identify the physical expectation value (A) in the infinite volume with lim (A; e, A) . 

L — >oo 

The solution of the implicit equation JSJ is obtained inverting the renormalization map 
R\ : e i— > e + K ( • ; e, A) , and corresponds to the sought-after interacting dispersion relation 

E(-,\,n) = R^(e-n)- (9) 
The expectation value (A) can finally be expressed as 

(A) x = lim E^(A;E(-,X,n),\') , (10) 

L^oo A'=A 

defined as a (formal) power series in A'. The power series expansion of (A}\ in the coupling 
constant A does not exist, since the power series expansion of 3^ L '(Jl, e, A') in the dispersion 
relation e does not exist. For A = 0, one recovers the expectation value of a system of free 
fermions. 

In the traditional Wilsonian approach to the renormalization group, the high energy de- 
grees of freedom are successively integrated out, modifying at each renormalization step the 
dispersion relation and the interaction. The bare Hamiltonian flows toward an effective Hamil- 
tonian, which should describe the low energy physics. Our construction begins directly with 
the interacting dispersion relation e, with the advantage of using the same propagator at each 
energy scale. The price to pay is the non-trivial inversion of the renormalization map. It is 
worth noting that the Hamiltonian {3} is not an effective Hamiltonian, but the starting point 
of the perturbation expansion. 

We now turn to the rigorous formulation of our results. Let 23 = {e : (p, fi) i— > e(p, fi)} be 
the space of all possible dispersion relations satisfying the following assumptions: (i) There is 
an open subset U of the first Brillouin zone, such that for all chemical potentials fi G /, where 
I is an interval of real numbers, e( • , jS) is smooth on U . (ii) For all fi £ I, the surface defined 
by is contained in U, and Ve(k) ^ for all k £ and [i e I. (iii) For all /iGl, the 
surface E 6iAI is strictly convex. We require that e(k) — /i satisfies the assumptions (i)-(iii). 

The first and second assumptions ensure the regularity and smoothness of the Fermi sur- 
face; in particular, van Hove singularities are excluded by (ii). The third assumption is used 
to control the value of Feynman graphs with two or more loops. The dispersion relations 
in 23 will correspond to the renormalized (or interacting) dispersion relations. Assumptions 
(i)-(iii) are hypotheses used in the proof of the following theorem on the existence of the 
thermodynamic limit: 

Theorem 1: Let e be a dispersion relation in 23, and v aT be a local interaction. Then 
there is a neighborhood of e in 23, and a renormalization map 



R[ L) : 



e i ► e + K^ L \--e,X) 



(11) 



such that the following holds. The expectation value of an operator A, defined by © is 
analytic in A, with a radius of analyticity that may shrink to zero as L — > oo. For L — > oo, 
the renormalization map tends to an invertible map R\ : e i— * R\e, defined as a formal power 
series in A, such that the limit 

E(A; e, A) = lim S (i) (A; e, A) (12) 

L— >oc 

is well-defined in the sense of formal power series, i.e., for ^ L \A\ e, A) = Xm>o C^S e )^ n > 

Ei L \A;e) L ^°°S n (A;e) (13) 
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for all n. In particular, the occupation number, defined as 

77 CT (k; e ,A) = lim S( L )(c+ a c kCT ;e,A) (14) 

L — >oo 

has a discontinuity as function of k on the Fermi surface E eA1 . The counter-term is determined 
by the condition that the self-energy vanishes order by order on the surface S e(J . 

A proof of this theorem will be given in a later paper [4] . The theorem 1 is a generalization 
for finite volumes of the main results of [5]- [8]. 

The proof of the analyticity of the expectation values in the finite volume would be trivial 
without the counter-term. One has to show that the discretized renormalization condition has 
a solution, which is analytic in A. Renormalization group ideas are used in order to control the 
value of Feynman graphs in the limit L — > oo. The vicinity of the Fermi surface Q is split into 
shells fitting into each other, the width of the shells getting smaller near the Fermi surface. The 
contributions of Feynman graphs at fixed scales are bounded by power counting, and careful 
summation over the scales shows the convergence of the Green's functions. When the width 
of the shell is of the order l/L, the power counting becomes worst, and scale decomposition 
proves to be useless. The regularized dispersion relation JJjJ provides an infrared cut-off at 
this critical scale. 

As a consequence of theorem 1, the density p(X, e) = lim p^ L \X, e), where 

L — >oc 

P (L) (A, e) = L- d ]T sW(c+ ffCkCT ; e, A), (15) 

k,(T 

is well-defined in the sense of formal power series. In the thermodynamic limit, the sum over 
quasi-momenta becomes on an integral over the first Brillouin zone, 

/d d k 
^p77 CT (k;e, A). (16) 

We prove the following statement: 

Theorem 2: The density l|16|l is independent of the interaction strength, 

p(A,e)=p(0,e) = 2Vol£ e , M , (17) 

denoting by Vol £ the volume enclosed in the surface E. 

Proof: First observe that, because of the regularized dispersion relation, the ground state 
is non-degenerate, providing the coupling constant is small enough. It follows that 
the expectation value of an operator A is given by its ground state expectation, S^ L ' (A, A) = 

(nW,io( £ )). 

Since the number operator N commutes with the Hamiltonian (J2J, is an eigenstate 
of both H ( L ) and N. The eigenvalues of N are positive integers, and the expectation value of 
the number operator is an integer. 

On the other hand, E( L \N, A) is analytic in A for small A. Hence 3^ L '(N,X) is constant 
with respect to A. In particular, p( L '(A, e) = S^ L '(A^, X)/L d is independent of A, and writing 
p^(X,e) — J2n>o Pn"\ e )^ n i we obtain p^{e) = for n > 1. Taking the limit L — > oo ends 
the proof of the first equality in Eq. (|17|) . 

The second equality in Eq. (|17fl is trivial: In a system of non-interacting fermions with 
a dispersion relation e(k, /i), the occupation number is 1, if e(k, p) < and vanishes for 
e(k, p) > 0. The second equality in Eq. ljT7l) follows then from Eq. (JTBJ) with A = 0. 
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We finally come to the proof of Luttinger's theorem. We now follow the original argumen- 
tation of Luttinger in [1] to conclude that, if the density is fixed to its physical value p p hys, 
the volume enclosed in the Fermi surface is independent of the interaction strength. 

In the free fermion approximation, obtained by setting A = 0, the (free) Fermi surface is 
given by 

So,^ = {k : e(k) - mo}, (18) 
where po is chosen such that the density of the physical system satisfies 

p phys = 2Vol S . (19) 

Consider now the system of interacting fermions. The occupation number of the system of 
many interacting fermions, given by n CT (k; p, X) = rj^Qs.; E{ ■ , A, p), A), has a jump on the 
surface . E( ■ , A, p) is the interacting dispersion relation defined in 10. 

By definition the physical Fermi surface is the surface on which the occupation number 
no-(k; fi, A) has a jump, and we can thus write 

S\, ^ = £ B ( . , m ,a)- (20) 

The density of the interacting system is given by p(X' , E{ • , A, p)), for A' = A. By the theo- 
rem 2, with . !fl> x) replaced by Sa, m and e( • , p) replaced by E( ■ , A, p), 

p(X', E(-,X, n)) = p(0, E(-,X, p)) - 2Vol S x ^. (21) 

In particular, the density of the interacting system is given by twice the volume enclosed in 
the physical Fermi surface. 

As in the case of free fermions, the chemical potential /i has to be adjusted to its physical 
value Pphys = p(X, pphys) in such a way that the density is fixed to its physical value, 

p(X, E(- ,A, Pphys)) = Pphys- (22) 

Equating the right hand side of Eq. lfT^|) with the left hand side of Eq. , and using Eq. l(2T|) 
leads to Luttinger's theorem, 

VolS^ phys =Vo\S 04l0 . (23) 

The interaction will in general deform the Fermi surface, and by the way change its volume. 
If the density is kept constant, and since the density is always given by the volume enclosed 
into the Fermi surface according to Eq. I|21|) . the chemical potential has to be modified in 
order to compensate for the change in the volume of the Fermi surface. 

As explained in the introduction, the present proof of Luttinger's theorem is perturbative, 
in the sense that only formal power series are considered. A non-perturbative proof of Lut- 
tinger's theorem, using the Fermi liquid construction of [3], is probably possible, under the 
hypotheses of [3]. However, it turns out to be very difficult to prove the invertibility of the 
renormalization map non-perturbatively. 

* * * 

A. Praz thanks C. Mudry for stimulating discussions, and careful reading of our manuscript. 
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